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ABSTRACT 

Context. The solar chromosphere and transition region are highly structured and complex regimes. A recent breakthrough has been 
the identification of dynamic fibrils observed in Ha as caused by field-aligned magnetoacoustic shocks. 
Aims. We seek to find whether such dynamic fibrils are also observed in Lya. 

Methods. We used a brief sequence of four high-resolution Lya images of the solar limb taken by the Very high Angular resolution 
ULtraviolet Telescope (VAULT), which displays many extending and retracting Lya jets. We measured their top trajectories and fitted 
parabolas to the 30 best-defined ones. 

Results. Most jet tops move supersonically. Half of them decelerate, sometimes superballistically, the others accelerate. This bifurca- 
tion may arise from incomplete sampling of recurrent jets. 

Conclusions. The similarities between dynamic Lya jets and Ha fibrils suggest that the magnetoacoustic shocks causing dynamic 
Ha fibrils also affect dynamic Lya jets. 

Key words. Sun: transition region - Sun: UV radiation 
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1. Introduction 

The solar chromosphere observed in Ho- is a bewildering mass 
of elongated features, but a breakthrough concerning so-called 
dynamic fib rils (henceforth H a DF s) has b ee n made in the 
studie s of De Pontieu et al. ('2004', '2007a'), Hansteen et alJ 
( |2006[) . jlouppe van der Voort et al. (2007), Heggland et al] 
( |2()07 ), and [L angange n et al.l (l2008alcl) following earlier work 
by ISuematsu et al., (,1995h . These studies have established that 
Ha DFs, which are rows of dark fibrilar features jutting out 
from plage and network with periodic extension and retraction, 
display repetitive mass loading by upward propagating magne- 
toacoustic shock waves driven by the global solar oscillations. 
Reduction of the effective gravity along tilted magnetic channels 
lowers their cutoff frequency and lets them propagate into the 
chromosphere, steepen into shocks, and repetitively lift the 
chromospheric-transition region interface. 

These studies are all b ased on optical observations, but 
Ide Wiin & De Pon tieu (2006) have studied transition-region jets 
in C IV in ultraviolet TRACE images and found remarkable mor- 
phological similarities between Ha DFs and their C iv counter- 
parts. The DFs may represent the injection of cool material pos- 
tulated by Judge (2008) as the source of hot sheaths making up 
the transition region. 

In this paper we report the presence of similar features in 
solar Lya images taken with the Very_high Angular resolution 
ULtraviolet Telescope (VAULT. iKorendvke et alj|2001h during 
its second rocket flight. VAULT acquired Lya images with much 
higher spatial and temporal resolution than the EUV imagery 
from the TRACE, SoHO, and STEREO sateUites but only for a 
fe w minutes. The Ly a fine structu re in these image s is dis cussed 
bv lPatsourakos etalj (|2007) and Judge & Centenol (|2008[) . Here 
we locate and study extending or retracting Lya brightness 



structures in images taken near the limb. We call these dynamic 
Ly a jets, abbreviated to Ly a DJs, and compare them to Ha DFs. 



2. Observations and measurements 

We used Lya imaged recorded during the second VAULT flight 
of 14 June 2002. The cadence was 17 s, the exposure time 1 s, 
the image scale 0.124arcsecpx"', and the angular resolution 
1/3 arcsec. The VAULT gratings isolated a 150 A wide spectral 
band around Lya and a narrow -band filter reduced the passband 
to 72 A FWHM. The resulting signal is about 95% pure Lya 
emission (Teriaca & Schiihle, personal c om munication). More 
detail is given in iKorendvke et alj ( l200lh and lPatsourakos et al] 
(I2007h . 

The VAULT-II flight recorded two Lya image sequences, a 
seventeen-image one stepping over an extended active area on 
the disk and a four-image one of quieter areas near the limb. We 
use only the four limb images because Lya jets are best observed 
as bright features in projection against the dark internetwork 
background. The four images were precisely co-aligned through 
cross-correlation using routines of P. Siitterlin. We then applied 
appropriate greyscaling to enhance their fine structure and used 
extensive visual inspection to identify Lya jets and to study 
their temporal behavior. They are similar to limb spicules in be- 
ing bright against dark internetwork and in showing hedgerow- 
clustering at network boundaries. Similar Lya jets are also seen 
in the disk sequence, but we decided to limit our measurements 
to the limb sequence that shows them best. It clearly shows tem- 
poral variation in the spatial extent of many Lya jets, even dur- 
ing the brief one-minute sequence duration. 



' The VAULT data archives are available 
I http : / /wwwsolar . nrl . navy . mil/rockets/vault/ | 



at 
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Fig. 1. Lya image taken at 18:16:56 UT covering the east limb. 
We invite the reader to inspect this image at large magnifica- 
tion with a pdf viewer. Lya jets are grouped in hedgerows at 
network boundaries that are best seen in projection against the 
dark internetwork. The rectangle selects the hedgerow enlarged 
in Fig. |2] The contrast is enhanced by logarithmic brightness 
scaling. Tickmark spacing: lOarcsec. 



The four images used here were taken between 18:16:56 
and 18:17:47 UT at the east limb. The first is shown in Fig. [1] 
Towards the limb it displays thick hedgerows of bright Lyo- jets, 
jutting out at network boundaries. They are remarkably similar, 
at reversed contrast, to the dark h edgerows in filtergrams t aken in 
the w ings of Ha (e.g.. Fig. 9.1 of lStixll2004l and Fig. 7 of Rutten' 
l2007h . in which Dopplershift of the line core into the wings se- 
lects dynamic fibrils over the quiescent network-spanning fibrils 
that dominate the scene at Ha line center. Such extended fib- 
rils appear dark in the VAULT-II images and provide the back- 
ground against which the hedgerows appear bright. Figure |2] 
shows an enlargement of one such hedgerow. The individual jets 
have widths of only one arcsec or less, requiring VAULT'S high 
resolution for identification and tracking. 

We have manually identified 55 Lya jets that exhibit mea- 
surable extension and/or retraction during the four-frame se- 
quence, all located in network clusters with convenient dark 
backgrounds. More weeding, described below, resulted in the se- 
lection of 30 of these Lya DJs for presentation. They are marked 




Fig. 2. Enlargement of the subfield outlined in Fig.[T] The rect- 
angle outlines the smaller subfield shown in Fig.|4]with Lya DJs 
3 and 1 1 . Tickmark spacing: 1 arcsec. 




Fig. 3. The upper part of Fig.[T]with 55 Lya DJs marked by plus 
signs. The 30 DJs selected for further analysis are labeled with 
their reference number Tickmark spacing: 10 arcsec. 



by their reference number in Fig. [3] Some examples are shown 
in Fig. m in the form of small cutout sequences from the four 
frames. There are many other features with similar morphology 
in Fig. |3j undoubtedly, many more might be recognized as Lya 
DJs in longer time sequences. However, there are also many jet- 
like features that do not change at all during the brief sequence 
duration. For example, the conspicuous jet to the right of retract- 
ing Lya DJ 24 in the bottom panels of Fig. |4]does not show any 
change. 

The cutouts in Fig. |4] illustrate that these Lya DJs appear 
rather fuzzy, making measurement of their extension or retrac- 
tion somewhat imprecise. At the suggestion of the referee to an 
earlier version of this paper, we developed a measuring method 
that includes error estimation. We first defined reference foot- 



J. Koza et al.: Dynamic Lyo' jets 




_^^^^^^^i^u 



Fig. 4. Time evolution of extending (nrs. 22, 19, 20, 3) and re- 
tracting (nrs. 11, 10, 9, 24) Lya DJs. The time increases with 
17-s intervals between the columns. The arrows define the jet 
axis. The tickmarks have 1-arcsec spacing. 



points and longitudinal axis orientations for each of the 55 Lya 
DJs to measure their lengths. Since many jets have no clear base 
due to crowding in the network from which they originate, the 
footpoint designation is rather arbitrary but this does not affect 
our length change measurements since we maintained the same 
footpoint for the three other frames after selecting it in the first 
frame. We similarly maintained the axis direction. No correction 
for projection effects was applied. The limbward viewing sug- 
gests that, on average, these Lya DJs are viewed more or less 
from aside. 

We then plotted the intensities along each jet against a length 
coordinate A measured in arcsec along the axis from the foot- 
point at A = 0. Figure |5] shows the intensity profiles for the same 
Ly Of DJs as displayed in Fig.|4] Each profile is box-car smoothed 
over 6 pixels (0.74 arcsec) to reduce noise. The differences be- 
tween the four curves in each panel show the temporal jet evo- 
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Fig. 5. Intensity profiles along the longitudinal axes of the Lya 
DJs shown in Fig. |4] The solid, dotted, dashed, and dashed- 
dotted curves correspond to the 1st, 2nd, 3rd, and 4th panel in 
the corresponding row of Fig. [H respectively, and show the tem- 
poral evolution. The thin-line markers specify the ranges of in- 
tensity thresholds defining each jet top. The reference footpoints 
with A = are fixed in time. 



lution. The DJs in the upper four panels extend; i.e., their outer 
parts become brighter Their lower parts sometimes weaken si- 
multaneously (as for DJ 22 in the first panel). The DJs in the 
lower four panels retract, without much change in their lower 
parts. Inspection of all 55 plots made us discard 15 DJs for not 
showing such regular progression in their outer parts. 

The question then was how to define the jet top locations in 
order to measure jet length variations in terms of their top dis- 
placements AAtop(f)- For each DJ we selected intensity thresh- 
olds at the onset of steep intensity decay with A and at the outer 
limit of visibility, as illustrated by horizontal dashes in Fig. |5] 
Picking an intensity value from such a range as the top loca- 
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Fig. 6. Trajectory measurements AAtop(f) of 30 jet tops along their axes. The heights of the symbols specify the error estimates. 
The solid curves are the parabolic fits. Nrs. 1-16: decelerating jets showing extension (1-8) and retraction (9-16). Nrs. 17-30: 
accelerating jets showing extension (17-22) and retraction (23-30). 



Table 1. Lya DJ excursion amplitudes and kinematic quantities. The first number of each entry is the median value, followed by 
the mean with standard deviation (between parentheses) an d the range (between angle brackets). The bottom row shows comparable 
determinations for Ha DFs from lDe Pontieu et al. I (l2007al) . 



feature type 


A [arcsec] 


\a\ [ms ^] 


v'max [kms"'] 


T [min] 


^ , decelerating 
Lya DJs 

accelerating 


0.6(0.8±0.5) (0.2,2.2) 


165 (235 ±173) (44,650) 
246 (257±175) (87,797) 


13 (17+9) (6,35) 
21 (18±7) (9,32) 


2.8 (3.3±2.2) (0.9,7.5) 
not defined 


Ho- DFs 




136 (146±56) (40,320) 


18 (18±6) (8,35) 


4.2 (4.8±1.4) (2, 10.8) 



tion in one profile yields samples of AAtop(0 for the other three 
by finding the displacement in A for the outermost pixel with 
that intensity. We automated this in a procedure using a random 
generator to pick threshold values from a normal distribution 
centered on and covering the selected range, and so obtained 
distributions of AAtop(f) displacements for each jet that yield 
rms error estimates. This was repeated using each of the four 
frames as initial reference; the four error estimates were aver- 
aged. These errors represent the threshold-definition uncertainty 
and the amount of gradient divergence between the four profiles 
per jet. They are smallest when the curves are parallel, as in the 
case of DJ 9 in the bottom-left panel of Fig.|5j 

Following th e exa mple of iHansteen et alj (l2006 l) and 

iDe P ontieu et^ ( 12007 a) we fitted the measured displacements, 
using the error estimates as weights, with parabolas AAtop(f) - 
vi (f - h) + (a/2) (f - fi)^ where vi is the onset velocity (positive 
when upward) at the time fi at which the first frame was taken 
and a the deceleration (when negative) or acceleration (when 
positive). The terminal velocity V4 at the time of the fourth frame 
is V4 = vi -I- fl (f4 - fi) where t^- t\ = 51 s. 

Since the onset and terminal velocities vi and V4 differ in all 
cases and can be negative, we define a maximum velocity v^ax 
as the higher of |vi| and |v4|. The apex of parabolic decelerated 
motion starting with v^ax is reached after v'max/|a| seconds. We 
estimate the total jet excursion duration T for decelerating jets 
by doubUng this value: T - 2vtnax/|a|- Since the four frames are 
likely to sample only part of a jet top trajectory, both v^ax and T 
represent only lower limits. 

Finally, we discarded 10 more DJs from our sample because 
the error estimates for their a and/or Vmax determinations ex- 
ceeded 100%. 



3. Results 

Figure|6]shows the top trajectory measurements, the correspond- 
ing error estimates, and the parabolic fits for the remaining 30 
Lya DJs. Table [T] specifies the average values of the trajectory 
amplitudes and the fit parameters together with their variances 
and ranges. 

Most jet tops do not travel far during the short sequence du- 
ration, only over a median distance of 0.6 arcsec, but these travel 
measurements are nevertheless significant thanks to VAULT'S 
0.3-arcsec resolution. Some of the decelerating jets have sig- 
nificantly greater deceleration -a than the solar surface gravity 
(274 ms~^); we call these superballistic. The maximum velocity 
estimates Vmax all but one exceed the chromospheric sound speed 
(8 kms"', e.g., Uitenbroek 2006). The excursion duration esti- 
mates T range from 1 to 7.5 minutes with a median at 2.8 min. 
Again, we emphasize that the last two parameters are likely to 
be underestimated when the brief image sequence did not cover 
the full jet-top excursion. 

Figure|2]plots the maximum velocity amplitude v^ax against 
the deceleration/acceleration amplitude \a\ for the selected 30 
Lya DJs. The scatter and the errors are large and the statistics 
low, but the plot suggests positive correlation between v^ax and 
\a\ for both the decelerating and the accelerating Lya DJs. 



4. Discussion 

Obviously, longer Lya image sequences at the VAULT-II resolu- 
tion are needed to gain better information than these four images 
provide. Nevertheless, they do show the existence of dynamic 
jets that fan out from network and extend or retract during the 
one-minute image sequence. 
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Fig. 7. Maximum velocity amplitude Vmax against decelera- 
tion/acceleration amplitude \a\ for the selected 30 Lya DJs. 
Filled symbols: decelerating DJs. Empty symbols: accelerating 
DJs. Dashed lines: chromospheric sound speed of 8kms ' and 
solar surface gravity. 



How do these Lya DJs compare to the now well-studied and 
fairly well-explained Ha DFs? They have similar shapes, occupy 
similar locations on the solar surface, show similar hedge-row 
morphology, and extend or retract along their length the same 
way. The measurements for the 16 decelerating ones in Table [T] 
show fairly good correspondence with the Ha DF measurements 
listed at the bottom of the table. Figure [T] suggests a positive cor- 
relation between v,nax and \a\ comparable to the linear relations 
for Ha DFs. Our estimate of the mean Lya DJ excursion dura- 
tion is about one and half minutes shorter than the average Ha 
DF periodicity reported by iDe Pontieu et al. I (I2007ah . but repre- 
sents only a lower limit. Three of the sixteen DJs show substan- 
tially greater superballistic deceleration than th e maxi mum of 
320 m s"2 for DFs reported by De Pontieu et a ll (l2007al) and the 
maxim um of 400 m s~^ repo rted by Roup pe van der Voort et alj 
(|2007|) and lDe Pontieu et al.1 (l2007b) for quiet-Sun mottles and 
for type-1 spicules, respectively, but with large uncertainty for 
one of them. We conclude that, overall, our decelerating Lya 
DJs are rather similar to Ha DFs. 

However, in contrast to the Ha DF studies, we found a 
roughly equal number of Lya DJs that significantly accelerate 
rather than decelerate. These are collected in the lower half of 
Fig. |6] The two sets of trajectories together suggest that they 
might actually be incomplete samples of sinusoidal motion or, as 
the referee has pointed out, sample the end of one shock-driven 
feature and the onset of the next one if these Lya DJs occur in 
succession as Ha DFs do. Note that our length estimates do not 
differ significantly between the two sets, but these are rather ar- 
bitrary. 

An obvious difference between Lya DJs and Ha DFs is that 
the former appear bright against dark backgrounds, the latter 
dark against bright backgrounds. Lya emissivity indicates the 
presence of gas that is sufficiently hot for collisional excitation 
of its upper level at 10.2 eV. Sufficient population of that level 
is also required to make f ibrils opaque in Ha. The recent study 
by iLeenaarts et al.l (l2007h suggests that cool fibrils may main- 
tain much greater Ha opacity than their temperature would sug- 
gest if they have recently undergone shock heating. Hence, Ha 
DFs may remain opaque also when they are cool during part 
of their lifetime. They will be much more opaque in Lya. Can 
the same fibril appear bright in Lya and dark in Ha? Perhaps 



through hot-sheath topology of the transition re gion, with the 
latter surrounding cool DFs as hot shells (e.g.. 'R uttenI 120071 : 
De Pontieu et al. 2007a). Th e neut r al-hyd rogen diffusion mech- 
anism recently proposed by Judge! (l2008l) seems a viable candi- 
date. Detailed time-dependent numerical MHD simulation may 
demonstrate this. However, realistic evaluation of Ha formation 
from such simulations remains a formidable challenge that may 
be helped by constraints derived from Lya. 

Finally, the Lya DJs studied here are likely related 
to chromosph eric phe nomena such as r egular spicules (e.g. 
ISterhng|l20cioh. straws (lRuttenll2006ll2007|). and tvpe-2 spicules 
jDe Pontieu et alJl2007bl;lLangangen et al.l l2 



l l2008bl) . 
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